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Abstract 

We consider three dimensional gravity with a positive cosmological constant and non- 
zero gravitational Chern-Simons term. This theory has inflating de Sitter solutions 
and local metric degrees of freedom. The Euclidean signature partition function of the 
theory is evaluated including both perturbative and non-perturbative corrections. The 
perturbative one-loop correction is computed using heat kernel techniques. The non- 
perturbative corrections come from gravitational instantons with non-trivial topology 
which can be enumerated explicitly. We compute the sum over an inflnite class of ge- 
ometries and show that, unlike the case of pure Einstein gravity, the partition function 
is finite. This demonstrates that the inclusion of non-trivial local degrees of freedom 
can render the sum over geometries convergent. 
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1 Introduction and Discussion 

Quantum gravity in de Sitter space is a notoriously difficult subject. Many basic questions 
- such as those involving the nature of observables in de Sitter space and the origin and 
interpretation of de Sitter entropy - remain unresolved. Many techniques which have proven 
useful in other circumstances, such as those involving AdS/CFT and supersymmetry, do not 
address these issues. It is natural therefore to focus on simple theories of de Sitter gravity 
where the deep questions of quantum cosmology can be discussed in a quantitative manner. 

Recently, we have focused on the question of three dimensional Einstein gravity with a 
positive cosmological constant [IJ. Although this theory contains inflating de Sitter solutions 
it has no local degrees of freedom. This allowed us to perform a series of precise computations 
which are impossible in more complicated theories of de Sitter gravity. We were able to 
compute the exact partition function in Euclidean signature, which is schematically written 
as a path integral of the form 

Z = [ Vge-'^^^ , (1.1) 
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where S[g] is the Einstein-Hilbert action. The physical interpretation of this path integral is 
the following. Euclidean path integrals are used to define states of the Lorentzian theory; the 
sum over compact geometries with specific boundary data defines the "wave function of the 



universe" in the sense of Hartle and Hawking [2] . The path integral (1.1) is an integral over 
compact manifolds without boundary and it is interpreted as the norm of the Hartle-Hawking 
wave function. 



Of course, functional integrals such as (1.1 ) cannot usually be defined precisely. However, 



one can give a precise definition to a path integral if one is able to 

• Ennumerate all solutions to the classical equations of motion, which appear as saddle 



point contributions to (1.1) 



• Compute the infinite series of quantum corrections around each saddle point 

• Perform the sum over saddle points, including this infinite series of quantum correc- 
tions. 

All three of these computations can be performed explicitly in the case of three dimensional 
Einstein gravity in de Sitter space. 

In fact, the first two of these tasks are not particularly difficult. The classification of 
solutions is related to the classification of spherical three-manifolds and proceeds much like 
the classification of crystallographic groups in three dimensions. The quantum corrections 
can be computed using the relationship between Einstein gravity and Chern-Simons theory. 
It is the third task - the sum over geometries - which proved most problematic. This sum is 
over an infinite number of topologically distinct geometries, and turns out to diverge in a way 
which cannot be cured using standard regularization techniques. Thus the Hartle-Hawking 
state of Einstein gravity is non-normalizable. 



In this paper we consider the path integral (1.1) for a slightly different theory of gravity, 
that of Einstein gravity with a positive cosmological constant and a gravitational Chern- 
Simons term. The classical solutions of the theory include the usual infiating de Sitter 
solutions of Einstein gravity with a positive cosmological consant. However, the Chern- 
Simons Lagrangian is third order in derivatives so the theory now possesses a local degree 
of freedom. This theory is known as Topologically Massive Gravity (TMG) and was first 
considered in [21 HI [5] . Despite the presence of this local degree of freedom, this partition 
function can be studied in considerable detail. In this paper we will focus primarily on those 
aspects of the computation which differ from the Einstein gravity case considered in fl^. 
Although this paper is self-contained we will quote certain results from [1]. 

The classical sum over geometries in TMG is described in section 2. The classification of 
Euclidean geometries is in fact identical to that of Einstein gravity. However, the classical 
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action of TMG evaluated on one of these solutions differs from that of Einstein gravity. 
The action is complex in Euclidean signature, so that each saddle point geometry is now 
weighted by a phase. If the Chern-Simons coupling is appropriately quantized these phases 
lead to cancelations in the sum over geometries. This has the effect of making the sum over 
geometries more convergent than in Einstein gravity. 

We then turn to the computation of quantum corrections, which are considered in section 
3. This is somewhat more complicated than in the Einstein gravity C3iSG, clS the theory 
possesses a local degree of freedom. Nevertheless, an explicit computation of the one loop 
determinant is possible. This is accomplished by using heat kernel techniques to compute the 
spectrum of the massive graviton wave operator and zeta function regularization to compute 
the regulated functional determinant. It also involves a careful treatment of the gauge fixing 
terms and Fadeev-Popov ghosts. In appendix A we give a careful derivation of the one-loop 
determinant of TMG using BRST quantization. In appendix B we give a detailed analysis, 
including a numerical study, of the resulting one loop determinant. 

Our conclusion is that, unlike the case of Einstein gravity, the sum over geometries 
converges when the gravitational Chern-Simons coupling takes certain discrete values. In 
this computation we include only those physically motivated saddles which have a natural 
interpretation in Lorentzian signature; this will be discussed in more detail in section 2. The 
inclusion of other saddles was considered in the case of Einstein gravity and shown not to 
quahtatively effect the result. We expect the same to be true here. We also include only the 
one-loop perturbative correction. Higher loop corrections are certainly present, and indeed 
quite difficult to compute in TMG. Our expectation is that these higher corrections will not 
effect our conclusions. This is based on our experience with Einstein gravity, where these 
higher order corrections were computed explicitly and shown not to alter the divergence 
structure of the sum over geometries [1]. 

We emphasize that our conclusion - that the sum over geometries converges only when a 
gravitational Chern-Simons term is included - is very similar to the corresponding result in 
Anti-de Sitter space. In that case the partition function of three dimensional Einstein gravity 
with a negative cosmological constant can be computed exactly, but the result does not have 
a quantum mechanical interpretation [6] . Once a gravitational Chern-Simons term is added, 
however, the sum over geometries has a natural quantum mechanical interpretation as the 
partition function of a dual conformal field theory [TJ |S] . We are finding a similar result in 
the case of a positive cosmological constant. This may indicate that pure quantum gravity 
in de Sitter space makes sense only when an appropriate gravitational Chern-Simons term 
is included. This is likely to have implications for the conjectured dS/CFT correspondence 
[0]; see e.g. [TUl |TT] for related considerations. We hope to return to this in the future. 

Finally, we note that other modifications of Einstein gravity to include additional degrees 
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of freedom may lead to similar results. Other straightforward extensions of three dimensional 
gravity include supersymmetric theories [I2l |13], generalized massive gravity [TH [15] and 
higher spin theories [161 EJ- It would be interesting to compute the partition function in 
these cases as well. 



2 The Classical Analysis 

In this section we study the partition function at the classical tree level approximation. In 



this limit the partition function (1.1) is given by its saddle point approximation 

Z = ^e-^°[^d+- . (2.1) 

9c 

Here the sum is over all classical solutions Qc to the Euclidean equations of motion and 
S'°[5fc] is the classical action. We will start by identifying the classical solutions Qc and 
describing their physical interpretation. We evaluate the tree level action 5° for TMG. We 
then explicitly perform the sum over geometries, including an infinite class of saddles with 
a clear physical interpretation in Lorentzian signature. 



2.1 Classical solutions 



Topologically Massive Gravity (TMG) is three dimensional general relativity with a gravi- 
tational Chern-Simons term [HIS]. Including a positive cosmological constant, the action in 
Lorentzian signature is 



5* 



IQtxG 



where Ics is the gravitational Chern-Simons term 



(2.2) 



(2.3) 



Here /i is a real coupling constant with dimensions of mass. The equations of motion of 
this theory are third order in derivatives of the metric, so unlike three dimensional Einstein 
gravity this theory possesses a propagating local degree of freedom. 
In Euclidean signature the action is 



S 



levrG 



/ (fx^(R 

Jm V 
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where Ics is now the Euclidean Chern-Simons term 

Ics 



(2.4) 



(2.5) 
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We emphasize the appearance of the factor of i in (2.4). This arises because the Chern- 
Simons Lagrangian is odd under time reversal t — )■ —t, so picks up a factor if i under the 
Wick rotation t — )■ it. Thus, as is usually the case for parity non-invariant theories, the 
action is complex in Euclidean signature. We will use units where i = 1 and define the 
dimensionless coupling k = i/AG. 



The equations of motion are found by varying (2.4) with respect to the metric. The 
Euclidean signature equations of motion are 

Gi,u + Qt^u = ifJ'~^Cf,u , (2.6) 

where C^i. is the Cotton tensor C^j, = In this paper we will restrict our attention 

to real solutions to the equations of motionj^ For a real metric, the left and right hand 



sides of (2.6) are purely real and purely imaginary, respectively. Thus they must vanish 
independently. In particular, the metric must obey the equation of motion of Einstein 
gravity with a positive cosmological constant 

G^. + 9t.u = . (2.7) 



When this equation is satisfied the right hand side of equation (2.6) will vanish automatically. 
We conclude that in Euchdean signature the equations of motion reduce to those of general 
relativity without a gravitational Chern-Simons termj^ 

It is now straightforward to enumerate the smooth solutions to the equations of motion. 
Equation (2.7) states that the metric must be locally S^. The classification of locally spherical 



geometries is a standard part of the classification of three-manifolds; see e.g. [19] for a review. 
We will simply summarize the results here. The solutions to the equations of motion are the 
three-manifolds S^/T, where F C 5*0(4) is a freely acting discrete subgroup of the isometry 
group of the sphere. There are an infinite number of possible subgroups F, which can be 
enumerated explicitly; they are central extensions of the crystallographic point groups in 
three dimensions. 

As discussed in pQ, there is a special class of solutions to the Euclidean equations of 
motion which have a natural physical interpretation. These are the lens spaces L{p, q), which 
are quotients of the three sphere by the cyclic group S^/Zp. These spaces are the positive 
cosmological constant analogue of the BTZ black hole solutions with negative cosmological 
constant [20]. They have a straightforward Lorentzian interpretation which we now review. 



^This represents a choice in our definition of the path integral as a sum over real metrics. This choice is 
justified by the fact that, as we will see later, the resulting partition function is convergent. However, other 
choices of integration contour through the space of metrics may be possible. This is the case in Chern-Simons 
gauge theory [18] so it would be reasonable to investigate a similar possibility here. 

^This argument was made in the case of TMG with a negative cosmological constant in [5]. 
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We start by considering the physics of a timehke observer in de Sitter space. This observer 
is in causal contact with the static patch of de Sitter space, which has metric 



ds^ = dr'^ — cos^ rdt'^ + sin^ rc 



(2i 



The Euchdean geometry is obtained by taking 



t ti 



it 



(2.9) 



which gives the metric 

ds^ = dr"^ + cos^ rrfti + sin^ rdcj)^ . (2.10) 

The geometry has to be smooth at r = 7r/2, implying that the Euclidean time coordinate 
must be periodically identified. Thus 



{tE, 4>) ~ (tE, 4>) + 27r(m, n) 



Vn, m G Z 



(2.11) 



The geometry (2.10) with the identifications (2.11) is the sphere S^. 

However, there are other identifications of the t and coordinates which make the ge- 



ometry (2.10) smooth. In particular 



+ 27r 



m q 

— , m — h n 

p p 



Vn, m G Z 



(2.12) 



is a smooth geometry provided p and q are relatively prime integers. These identifications 
define the lens space L{p, q) = S^/Zp] the sphere is L(l, 0). We note that a shift of q by 



a multiple of p can be absorbed into a shift of n in (2.12). Therefore the parameter q can be 
taken to be between 1 and p. 

To understand the physics of these geometries, we note that the operator that generates 



the identification (2.12) is 



p = exp H iqj] , (2.13) 

\ P V J 

where the charges H and J generate time translation and rotation, respectively. Equation 



(2.13) can be regarded as a density matrix which defines a grand canonical ensemble with 



temperature (3 = 2Tr/p and angular potential 9 = 2Tiiq/p. 

This provides a natural physical interpretation for the lens space L{p,q). At the level 
of quantum field theory in a fixed de Sitter background, correlation functions in L{p, q) 
can be Wick rotated to obtain correlation functions in de Sitter which are evaluated in a 
thermal state of fixed temperature and angular potential. Once gravitational effects are 
included, each L{p, q) represents a contribution to the Hartle-Hawking state. The dominant 
contribution is from the sphere = L(1,0), which has vanishing angular potential. This 
contribution describes the standard thermal behaviour of de Sitter space. The other lens 
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spaces give subleading contributions which lead to deviations from thermahty. The partition 
function is then the norm of this Hartle-Hawking state. 

We emphasize that the lens spaces are the only Euclidean geometries with a clear 
Lorentzian interpretation. The other saddles of the form S^/T, where F is not a cyclic 
group, can not be Wick rotated to the static patch of de Sitter space. We will therefore 
focus in what follows on the sum over lens spaces. The sum over other saddles can be in- 
cluded, but we do not expect that this will lead to qualitatively different results. This was 
discussed in detail in [1]. 



2.2 The Gravitational Chern-Simons Action 

We now need to compute the classical action of one of the L{p, q) saddles, including the effect 
of the gravitational Chern-Simons term. This is most easily done using the Chern-Simons 
formulation of three dimensional gravity, where the action is expressed not as a function of 
the metric but rather as a function of the frame fields e"^ and the connection u"". We will 
follow the same conventions as [I]. 

We first note that the Chern-Simons formulation of TMG is somewhat more subtle than 
that of Einstein gravity. In the Einstein gravity case, the equations of motion are completely 
equivalent to those of a Chern-Simons theory; this is implied by the famous equivalence 
between the second order (metric) and first order (Palatini) formulations of general relativity. 
In the case of TMG, however, the equations of motion are third order in the metric and 
describe a propagating local degree of freedom. This local degree of freedom is not present 
in the Chern-Simons formulation; Chern-Simons theory is topological. So the two theories 
are not equivalent, even classically. 

Nevertheless, the Chern-Simons formulation can be used to evaluate the classical action of 
TMG for certain solutions. In particular, if the Cotton tensor vanishes (i.e. we are studying 
a solution of Einstein gravity) then the TMG action is precisely that of a Chern-Simons 
gauge theory. This is exactly what happens when we restrict our attention to real metrics 
in Euclidean signature. 

Explicitly, the Chern-Simons action is 

I[A] = j 1i(^AAdA+'^AAAAA^ , (2.14) 
where A is an SU{2) connection and Tr is the usual trace on the SU{2) Lie algebra. When 



the Cotton tensor vanishes, the TMG action (2.4) can be written as that of an SU{2) x SU{2) 
Chern-Simons theory 

S = - (l^/Kl + ^ . (2.15) 
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The gauge fields A± are related to the frame fields and spin connection cu" = ^e°'fy'^u'^c by 



Al = Lo''± e" 



A± = AlTa . 



(2.16) 



Here are the SU (2) generators. The levels k± are complex and are related to the gravi- 
tational couplings by 



i 

2G 



2k 



(2.17) 



We now use this to compute the action of a lens space. 

Using the metric (2.10), the SU{2) x SU{2) connection on L{p,q) is 

A± = A%Ta = ±Ti dr + (cos rTa ± smrT3)de± , 

where 9± = (p ±tE obey the identifications 

n(q ± 1) — mp , , _ 
e±^e± + 27T^ ^ , Vn, m e Z . 

P 



(2.18) 



(2.19) 



We note that the lens space L{p,q) = S^/Zp contains a non-contractible cycle given by 
the Zp quotient; this is the cycle {n,m) = (1,0) in (2.19). To describe this connection 
more geometrically, we should compute its holonomy around this non-contractible cycle. As 
7ri(L(j9, g)) = Zp, this holonomy is a p*'* root of unity in SU{2) x SU{2) so is conjugate to 
an element 

2iTin+/p n \ / 27rin_/p f) \ \ 

e-wJ' e--«-/J^^^(')"^^(') ('-'^^ 



It is straightforward to compute the holonomy of the connection (2.18) and show that it is 
given by 

' q + l q-l 



in+,n_) 



(2.21) 



2 ' 2 

We now compute the action of TMG using known expressions for the SU{2) Chern- 
Simons invariant on a Lens space. For an SU (2) connection on a lens space with holonomy 
n, the Chern-Simons invariant is 



1 

87r2 



Ti { AAdA + -AAAAA 



T 2 

— n 
P 



where q* is the inverse of q mod p: 



q*q = 1 modp 



(2.22) 



(2.23) 



From ( |2.21[ ) and ( l2.22[ ) we find that the action of TMG is 

'ik. 



(P>9) 



exp , —I[A^] + '^I[A_ 



q* q ((f' + 1)5'* 

exp ( 7rz(A;_|_ — k ) h 7Ti{k+ + k. ^ 



p 2p 

[2'Kk .k{q + q*)\ ,^ 
= exp ^m-— — — . (2.24) 

VP P J 

The first term in this action is real; it is simply the usual Einstein action. This term is 
proportional to the volume of the lens space L{p,q) = S^/Zp. The second term is purely 
imaginary and comes from the gravitational Chern-Simons term. 

We note that in general the action of a Chern-Simons gauge theory is invariant under 
large gauge transformations only if the real part of the Chern-Simons coupling (i.e. the level) 
is an integer. In the present case this implies that the gravitational Chern-Simons coupling 
must be quantizec^ 

A e Z . (2.25) 

In the gravitational language this condition is necessary for the action of TMG to be invariant 
under large diffeomorphisms in Euclidean signature. From now on we will demand that /i is 



quantized in accordance with (2.25). Indeed, it is only with this quantization that the action 



(2.24) is invariant under shifts of q or q* by a multiple of p. 



2.3 The Sum over Geometries 

We now compute the sum over geometries including the tree level contributions described 
above. 



The sum over lens spaces is 



^lens 



p=i (p,g)=i 



= ^e2-'=/^5(A;/(2/x),A:/(2/.),p) . (2.26) 
p=i 

Here S{a, b, m) is the Kloosterman sum 

m 

S{a,b,m) = exp (27ri (an* + 6n)/m) , (2.27) 

n = l 

(n ,m) — l 

where n* is the inverse of n mod m. Expanding the exponential we find 



00 



Ziens = I J2p-''S{k/{2f,),k/{2fi),p) ] . (2.28) 



T ! 

r=0 Vp=l 



■^We note that a factor of two appears in this expression because the Chern-Simons gauge group corre- 
sponding to three dimensional gravity not actually SU{2) x SU{2) but rather 50(4) = SU{2) x SU{2)/12. 
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The quantity in parenthesis defines what is known as the Kloosterman zeta function, which 
appears frequently in number theory. We include in appendix [C] a summary of various 
features of this zeta function. This zeta function has a pole at r = 0, implying that the sum 
over tree level geometries diverges. 

We note that the divergence described above is not of particular physical significance, 



as one-loop effects can (and will) change the p dependence of the terms in the sum (2.28). 



However, it is useful at this point to note that the the sum (2.28) is already considerably 
more convergent than the corresponding sum in Einstein gravity. In that case all of the r < 2 
terms diverged |T]. The improved convergence comes from the phases which appeared in the 



action (2.24), which lead to cancelations in the sum over geometries. These cancelations 
appear because the phase of the Kloosterman sum is essentially randomly distributed]^ 
We now turn to the discussion of quantum corrections, which will render this sum finite. 

3 Quantum Corrections 

In this section we compute the one-loop quantum corrections to the TMG partition function 
by evaluating the appropriate functional determinants. We will then perform the sum over 
geometries including quantum corrections and demonstrate that the result is finite. Technical 
details related to the computation of one-loop determinants are included in the appendices. 

3.1 One loop determinant 

The one loop contribution to the partition function 



Z= / _^e-^[^] = ye-^°+^^+-, (3.1) 
J Vdis 

9c 

is obtained by integrating over the linearized fiuctuations h around each classical saddle. 
The computation is rather similar to that in Einstein gravity, the only difference being that 
we now have an additional propagating degree of freedom. 

A detailed derivation of the functional determinants appearing in 5*^ is provided in ap- 
pendix |A| This involves a careful accounting of the residual gauge symmetry and Fadeev- 
Popov ghosts in the measure factor. Related discussions in the case of Einstein gravity 
appear in [211 ISS |26] and for topologically massive gravity with a negative cosmological 
constant in [27]. 

The answer is the following product of functional determinants 

m det'^/2(-A(i) -2)r , , 

det'i/2(l + z/i-iDM)TT det'i/2(-A(2) +2)tt 



^In fact, the distribution of the values of the Kloosterman sum is a topic of great interest to number 
theorists [22ll23l. 
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Here A(j) = V"Vq, is the Laplacian acting on a field of spin j and tlie subscript T (TT) 
refers to the functional determinant on the space of transverse vectors (transverse traceless 
tensors). The prime indicates that we restrict to the positive part of the spectrum and Dzm 
denotes the contribution from zero modes. The operator Dm acts on symmetric 2-tensors 
and is defined by 

DmT^, = ^(e/^V„T,^ + eJ'^VaT^p) . (3.3) 

As one might expect, this formula is closely related to that of Einstein gravity. The only 
difference is the determinant involving Dm, which comes from the local degree of freedom. 



Thus we can write (3.2) as 



7(1) 7(1) 



(3.4) 



where Z^^^ is the one-loop determinant for Einstein gravity and 



.(1) 



'MG ~ V-^ "T ''P J^M)TT ■ (3.5) 

We may now use the results of [1], where Z^)^ was evaluated explicitly. For a lens space 
L(p, q) with g 7^ ±1 mod p we have 

cos I — COS COS 



7(1) 

^Ein.lens 



27r 



P 



and 



'(1) 



P 



7(1) 

Ein,(p,p-1) 



TT 



P ) 



271 



COS 



27rg* 



P 



Ein,(p,l) 

for p ^ 2. Finally, for and L(2, 1) we have 



2kp' 



sm 



P 



7(1) 

^Ein.S^ 



¥k 



7(1) 

^Ein,(2,l) 



TT 



2^^k 



(3.6) 



(3.7) 



(3.^ 



We now turn to the computation of Z^^^q. 



3.2 Massive graviton determinant 

We need to compute the one loop determinant for the massive mode 



7(1) 



det-^^\l + ifi-^D 



TT 



(3.9) 



We will use heat kernel techniques to compute the spectrum of this operator, and zeta 
function regularization to compute the functional determinant. This regularization proceeds 
as follows. The operator Dm has discrete eigenvalues A„ with degeneracies dn which we will 
compute explicitly below. The logarithm of the functional determinant is 



^^rf„ln(l ^A„) . 



(3.10) 
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This can be written in terms of the zeta function 



Cmg{s) = 5Z n 



dr. 



as 



log 



1 d 



2ds 



(3.11) 



(3.12) 



The functional determinant is formally divergent, but by analytically continuing Cmg{s) to 
the whole complex s plane and evaluating (3.12) at s = we obtain a regulated answer for 
the determinant. 



We now need to construct the zeta function (3.11) by identifying the eigenvalues and 
degeneracies of Dm- The first observation is that the square of Dm is the standard laplacian 
when acting on symmetric, transverse traceless tensors 



n2 ^ (TT) 



A(2)+3)T 



(TT) 
fit/ 



(3.13) 



This follows directly from (3.3). This identity relates the spectrum of A(2) with the spec- 
trum of Dm- If the operator —A (2) has eigenvalues a„ with degeneracies dn, then Dm has 
eigenvalues A„ = A/a„ + 3 with the same degeneracy dn^ 

Using this relation between Dm and A(2) we can build Cmg{s)- From [28] we know 
the spectrum of the spin-2 Laplacian on L{p,q). The eigenvalues and degeneracies of the 
transverse-traceless modes are 



and 



Here n 



and 



dr, 



0,1, 



mGZp 

and we have defined 
27r 



r\ - T2 



P 



Ti + r2 



27r 
P 



+ 



Xn{r) 



sin((2n + l)f: 



sm 



Gathering the above results the zeta function (3.11) becomes 

00 



(3.14) 
(3.15) 

(3.16) 
(3.17) 

(3.18) 



n=0 



^In principle there is an ambiguity in the sign of A„, since we could take either branch of the square root. 
Because Dm is a self-adjoint operator its spectrum is bounded, so the sign of A„ is fixed for all but a finite 
number of eigenvalues. This amounts to an ambiguity of the phase of the determinant det(l + i^^^Dm)^^'^ 
coming from this finite product of eigenvalues. At the end of the day we will only be interested in the norm 
of the determinant, so this ambiguity is irrelevant for our purposes. 
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We now wish to understand the analytic properties of (3.18). Even though it is difficult 
to evaluate explicitly equation (3.15) for the coefficients dn we note that the sum can be re- 
arranged in terms of Hurwitz zeta functions. This follows from the fact that the coefficients 
dn are almost periodic 



d. 



rp+J 



dj + 2r{rp + 2j + 6) 



(3.19) 



for g 7^ ±1 mod p. Special values need to be discussed separately, which we will do in 



appendix [Bj This allows us to write (|3.18|) as 

Cmg(s) = 



/i 

ip 

+2p 



[dj -2p\aj\'^]C{s,aj) 



(-) E 



C(s - 2, aj) + -ifiCis - 1, aj) 



where 



p 



(3.20) 



(3.21) 



The advantage of this expression is that we know the analytic properties of the Hurwitz 
function C{s,a) and its derivatives with respect to s. Using (C.5) and (3.20) we write (3.12) 

as 



^l^i) + ^ E - 2pl«^ n [l^r (a,) - I ln(27r)] 



C'(-2,aj) + -i/iC'(-l,aj) 



(3.22) 



It is also useful to compute the norm of Z^^^q, which 'n 

2Re{N) p-1 



7(1) 12 



(2vr) 



-A/2g7rIm(Af)+B /^/^A 



Pj 



n 

j=0 



r 



j + 3-ifi 
p 



(dj-2p|a,|2) 



(3.23) 



In (3.22) and (3.23) we have introduced 



A = J2 {d, - 2p\a,\') 



B 



Re('{—2,aj) film('{—l,aj] 



(3.24) 



^This is a variant of the Chowla-Selberg formula. If (dj — 2p|ajp) is of order one we can easily compute 
the product of gamma functions using Gauss's multiplication formula. 
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3.3 The Sum over Geometries 



We are now in a position to compute the sum over geometries, including the quantum 
correction computed above. The partition function is 



E 2 

(P,9) = l 



(0) ry{l) 



(3.25) 



where z'^^^-^ is given by (2.24) and the one-loop corrections are 



.(1) 



7(1) 7(1) 



(3.26) 



which are given in (3.6), (3.7), (3.8) and (3.22) 



Our claim is that (3.25) is finite, unlike the Einstein gravity case. To see this, we first 
investigate the sum excluding the one- loop contribution of the massive graviton z'^^^q. The 
sum includes various terms (such as those coming from Lens spaces with g = ±1 mod p) 



which are finite when summed over p. The only term which is not finite comes from (3.6) 
and is 



E 

{p,'?)=i 



7(0) 7(1) 

^(p,g)^Ein 



27r 



E^^''^'^'W(2/i),fc/(2/i),p) + .. 



p=i 



k~'Y.^^[Y^p-^'^'^S{k/{2^,\k|{2^,U)] +... (3.27) 



r=0 



Again we have encountered the Kloosterman zeta function, which is the quantity in paren- 
thesis in the second line. As described in appendix C, this zeta function is finite for every 
value of r > 0. The only potentially problematic term is the one with r = 0, where we are 
considering the sum 



X^j9-2^5(A;/(2/i),fc/(2/x),p) 



(3.28) 



p=i 



with s = 1/2. This particular value of the Kloosterman zeta function is of considerable 
mathematical interest, as the behaviour of the Kloosterman zeta function on the line = 
1/2 is related to the spectrum of the hyperbolic laplacian on the modular surface 'H/SL{2, Z). 
In general, the behaviour of sum (3.28) at s = 1/2 is a difficult number theoretic question; 
we refer the reader to [29l [231 for more details. 



In the above discussion, we have neglected the important Z^^^^ term. To include the 



effects of this term we must understand its p dependence, which is not immediately evident 
from (3.22). In order to demonstrate that the partition function is finite it is sufficient 



to argue that for large values of p the determinant decreases sufficiently quickly as a 
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function of p. This can be checked numerically, as we describe in appendix |B| | In particular 
we note that for every value of ^ there is an e > such that 

2Re(7V) p-1 

n 

j=Q 



7(1) |2 



< e 



Trlm{N)+B 



Pj 



r 



J + 3 
P 



(d,-2p|a,f) 



o{p- 



(3.29) 



for sufficiently large p. Even without analytic expression for z9:)q as a function of p and q 



the upper bound (3.29) is sufficient for our purposes. Armed with this result we can now 



reconsider the troublesome r = term in equation (3.27). In this case we are now considering 



a zeta function of the form (3.28) at s = 1/2 + e for some positive value of e. This sum is 



convergent. We conclude that the partition function of topologically massive gravity in de 
Sitter space is finite. 
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A Perturbative analysis of TMG 

In this appendix we construct explicitly the one-loop determinant of TMG. We start by 
expanding the metric around a background solution g of the equations of motion 

g^lu = g^iu + hf,^ , (A.l) 

where /i is a linearized metric fluctuation. The action takes the form 

S[g + h] = S'\g] + ^(^[/i] + S„[h] + Sghost + Oih") . (A.2) 

Here is the on-shell action, S^^^k is the action (2.4) expanded to quadratic order in h and 
'S'gauge and S'ghost are the contributions from the gauge flxing terms and ghost determinant, 
which will be computed explicitly in the next subsection using BRST quantization. 

To obtain St^^^^ it is convenient to decompose the metric perturbation into its scalar trace 
and traceless tensor parts 

(A.3) 



with 



0|^ The curvature tensor of the background metric satisfles 



R 



6 



R = 6 



(A.4) 



^Although we give only numeric evidence in appendix B, we suspect that an analytic proof is possible. 
®In this section indices are raised and lowered with g^^ and all covariant derivatives are with respect to 
the background metric. 
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so the quadratic action can be written as 



9 

(A.5) 



Here Aq) = V"Vo is the Laplacian acting on a field of spin j and we have introduced the 
operator Dm 



1 
2 

which acts on symmetric 2-tensors. We have defined the fields 



DmT^. = 7^{e^^No.T,p + e,"^V,T^^) , (A.6) 



V = (A.7) 

A.l BRST quantization 

We now turn to the gauge fixing conditions and Fadeev-Popov determinants. For a higher 
derivative theory such as TMG it is convenient to use the BRST approach, which gives a 
derivation of the ghost determinants while preserving gauge invariance. For a more detailed 
discussion of the BRST quantization of Chern-Simons theories and TMG see respectively 
[30] and [3l]. 

The BRST transformations are 



5c" = -B" , 

SB" = , (A.8) 

where c and c are the Fadeev-Popov ghost and anti-ghost vector fields and B is an auxiliary 
field. The notation is chosen to resemble the usual Fadeev-Popov procedure, but we note 
that c and c are not complex conjugates. It is straightforward to check that 6^ = 0. The 
gauge fixed action is 

Stot = S + 6<i/, (A.9) 

where \E' is an arbitrary functional of ghost number -1. The variation 6"^ includes both S'gauge 
and Sghost- The traditional choice for \l/ is 

^ = / d'x^c'' (\b^ - G,] , (A.IO) 
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where Ga is the gauge fixing condition. One can then solve for B and plug back into the 
action. However, this procedure gives us a gauge-fixing term proportional to GaG°' whereas 
we need a term third order in derivatives. So we will take the functional \I' to be 



where 



^ I d'x^-c"iDs)£(lB,-G, 



/3 ' 



(A.ll) 



(A.12) 



Since we are interested in the one-loop part of the path-integral, we expand the metric around 
a saddle point 

9tiu = 9tiu + h^y , (A. 13) 

and keep the terms that are at most quadratic in 5, c, c and h^^. We will take the gauge- 
fixing condition to be 

G^ = V^{Kp-]^g^ph) , (A.14) 
where h = g°'^hai3. At linear order 

SG^^^ = (g^f^V^ + R^^y , (A.15) 



so 



(A.16) 



Rewriting the equation above gives 

S^^'^ = ^ld'xV^(^-^-G''{DB)apG^ + (^{DB)^'ig,pV' + R,py^ 

+ 16^ / '^'"'^ l^^" " G^iDsUiB^ - G^) . (A.17) 



Using (|A.14|) and (|A.1|) we obtain the following form for the gauge- fixing term 

1 



32nG 
1 



V + 2c"(Db) J(^77/3V' + Rj^y + - G'°)(DB)a/3(5^ - G^) ) (A.18 



A. 2 One-loop Partition Function 



Collecting all the terms quadratic in h from (A.4), (A.5) and (A.18), the one- loop effective 
action is 
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1 



327rG 
1 

'l67rG 



(1 + ^f^-'DM) (-A(2) + 2)0^, + ^0(A(o) + 4)0 
-c-{Ds)f{-A^i) - 2y - \{B- - G-){DsU{B^ - G^) 



(A. 19) 



The one-loop partition function is 



V(f)f,^V(f)VcVcVB e 



-50+51 



det(DB(-A(i)-2))det-i/2p^) 



det^/2 + i^^-Wm) (-A(2) + 2)) deti/2(-A(o) - 4) ' 



(A.20) 



One can further simphfy this by decomposing these operators into transverse modes: 



det((l + ^/x-iDM)(-A(2) + 2)) 
det(DB(-A(i) - 2)) 



det((l + z^-iDm)(-A(2) + 2))t det(D5(-A(i) - 2)) 
det(DB(-A(i) - 2))t det(-A(o) - 4) . (A.21) 



Then, 



det^/2 ((1 + ifi-WM) (-A(2) + 2))^^ 



1/2/ 



-'Dm) (-A 

where the subscripts denote that the operators act on transverse and traceless tensors 



(A.22) 



We now proceed by sphtting the determinants appearing in (A.22) using the formula 
det(74i?) = det(y4) det(-B). While this splitting is straightforward for finite dimensional op- 
erators, a subtlety arises for infinite dimensional operators of the sort appearing here. The 
problem is that while the expressions det(y4i?) and det(yl) det(-B) are equal when regarded 
as formal products of eigenvalues, the zeta function regularizations of these infinite prod- 
ucts may differ. The difference between these two expressions is known as a multiplicative 



anomaly; see e.g. [321 [33] for more details. The operators of interest in (A.22) are all powers 
of the Laplace operator. The multiplicative anomaly for Laplace operators has been shown 
to vanish in odd dimensions [32] . We will therefore proceed to split the determinants without 
including a multiplicative anomaly, although this issue may be worthy of further study. 



Finally, we note that some of the operators appearing in (A.22) may have zero or negative 



eigenvalues. These arise when one studies geometries which possess isometries or conformal 
killing vectors. To treat the zero modes properly one should not include them in the deter- 



minant (A.22), but instead integrate directly over the appropriate moduli space of collective 
coordinates. The negative modes (which come from the scalar operator A(o)) are treated by 
rotating the contour of integration in field-space. Both of these subtleties appear in Einstein 
gravity and were studied in detail in pQ; the analysis in TMG is identical, so we will not 
repeat it here. 
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We are left with our final expression for the one-loop determinant 

det'V2(_A(i) - 2)t 



(A.23) 



det''/'(l + ifx-WM)TT det''/2(-A(2) + 2)tt 
Here Dzm is the zero-mode piece. The prime indicates that we restrict to positive eigenvalues. 



B More on the massive graviton 



In this appendix we discuss some of the properties of the determinant (3.5). We will start by 



evaluating Z^^^^ for special values of (p, q). In the next subsection we bound the determinant 
using numerical estimates. 



B.l Special cases of L(p,q) 



The zeta function associated to Z^^j^q is 



(B.i: 



n=0 



with dn given by (3.15). 



We start by considering the three-sphere 5*^ = L(l, 0). The degeneracies for this case are 



dn = 2{{n + 3Y - 4) 



(B.2) 



This allows us to write (B.l) as 



Cmg,53 = ^i-ifiY [C{s - 2, + 2zK(s - 1, -^/i) - + 4)C(s, -z/i)] + — 



+ z/i" 



[B.3) 

where we neglect terms independent of s. From here it is straightforward to compute Z^^^q 



using (3.12). The resulting expression is not illuminating, hence we leave it as an exercise 



for the curious reader. 

The next geometry that needs separate treatment is L(2, 1). Here, the degeneracies are 



c?2„ = 2(2n + l)(2n + 5) , d2n+i = 



(B.4) 



and the zeta function is 

Cj\/G,(2,1) = < 



B.5) 
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B.2 Numerical estimates of Z 



(1) 

MG 



Here we numerically explore the p-dependence of the expression (3.23) for the determinant 
z'^Dq. The claim is that for large values of p - in particular p ^ jj,- this determinant decreases 
as a negative power of p. Our starting point is the inequality in (3.29), which reads 

2Re(7V) p-1 



7(1) |2 
'MG\ 



< e 



7rIm(7V)+_B 



pj 



n 

j=0 



r 



J + 3 
p 



(B.6) 



The coefficients aj, N and B are given in (3.21) and (3.24). 
We first focus on the exponent in the gamma functions. 



p-i 

n 

j=0 



r 



J + 3 
p 



{dj-2p\a,\^) 



(B.7) 



we 



The coefficients dj for L{p,q) are bounded by the values (B.2), and from figure B.2 
observe that the maximum value of dj is of order p which occurs at j = p — 1. For sufficiently 
large values of p, {dj — 2p\aj\'^) is always negative which makes the gamma function in (B.6) 
negligible. 

Next, we turn to which is defined in ( |3.24 ). The real part of is 



ReN 



E 



<'i + + 3)' + 1)0 + 3) 



3p 



6 p 



which is strictly positive, whereas the imaginary part is 

5 (2/i2 + 37) 
2 6p 



ImiV = /i(^l5^rf, 



p 
3 



< 



2 



Therefore the asymptotic behavior of terms in (B.6) that depend on is 

2Re{N) / N 2|Re(Af)| 



,7rIm(Af) 



p) 



~ e 



P 



+ . . . 



(B.8) 



(B.9) 



(B.IO) 



The final term we need to consider is B which is considerably more difficult. Taking 
p^ li^Q have 



B 



4P E ^' -2' 



p 



:C'(-2) 



0.12 

p 



(B.ll) 



Hence < 1 and can be neglected. 

Gathering our estimates and observations above, we conclude that for every value of fi 
there is a positive number e such that 



1^(1) |2 < p7TtM/2 



p) 



2Kc{N) p-1 

n 

j=0 



r 



J + 3 
p 



(dj-2p|a,|2) 



0{p' 



(B.12) 



We have extensively checked this bound numerically. While it is possible to make much 
stronger estimates for Zmg, equation (B.12) will be sufficient for our purposes. 
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Max djip,q) 



200 ^ 
100 - 



Figure 1: The plot shows the degeneracy of the highest degenerate state as a function of p. 
It suggests a hnear growth in p with different slopes for even and odd p. 



C Dirichlet Series 

In this appendix we summarize some useful number-theoretic formulae. 

Riemann zeta function: 

The Riemann zeta function ({s) is the analytic continuation of the series 

c(^) = E^= n i^-pr\ (c.i) 

n=l p prime 

to the complex s plane. The function has a simple pole at s = 1 and Laurent series 

k=0 

where jk is the Stieltjes constant. Some useful values of ({s) are 

C(0) = -^, ^C(0) = -^ln(27r). (C.3) 

Hurwitz zeta function: 

A simple generalization of the Riemann zeta function is the Hurwitz function 

oo ^ 

^ in + a)'' 

It is a meromorphic function in s and 5R(a) > —1 with a simple pole at s = 1. We will need 
the following values 

C(0, a) = ^ - a , ^C(0, a) = \n{T{a)) - \ ln(27r) . (C.5) 
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so that in particular 

^C(0, a) + ^C(0, -a) = - ln(sin(7ra)) - ln(-2a) . (C.6) 
as as 

Kloosterman Zeta Function 

We now summarize a few features of Kloosterman zeta functions. The Kloosterman sum 
is defined as 

m 

S{a,b,m) = exp{27ii{an* + bn) / m) , (C.7) 

n=l 
(Ti,m) — 1 

where n* is the inverse of n modulo m. We are interested in sums of the form 

oo 

L{m,n;s) = ''^^p^'^^S{m,n;p) . (C.8) 

p=i 

This is known as the Kloosterman zeta function (see [29J for details). This series converges 
absolutely when 3fJs > 1/2. 

The analytic properties are most conveniently summarized by the function 

Z{m,n;s) = y^p~^S{m,n;p)J2s-i{ — Vrrm) , (C.9) 

p=i 

when mn positive, with a similar formula for mn negative. Using the Neumann expansion 

A:=0 

with z = —Jmn and z/ = 2s — 1 we see that 

p V 

22^v^ (2(s + A;) - l)r(2s - 1 + A;) ^ , ,^ ^ , 

L m,n;s = ^^.-i >^ — ' fp -Z{m,n-s + k) . C.ll 

The only poles of L{m,n] s) on the real s axis come from the gamma function, which 
has simple poles at the non-positive integers. For the k = term in the sum these poles are 
cancelled by the coefficient 2{s + k) — 1. Thus L{m, n; s) has no pole at s = 1/2. However, 
when s = —n/2, n = 0, 1, . . . there will be simple poles. For example, there is a pole at 
s = with non-zero residue coming from the k = 1 term: 

L(m,n] s) ^ -AnmnZ {171,^,1) + .. . . (C.12) 

s 

Similar conclusions hold for the case where mn is negative. 
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